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I. INTRODUCTION 

Quantum non-demolition (QND) measurements enable 
repeated measurements on quantum objects with accu- 
racy levels exceeding the standard quantum limiti. Such 
QND measurements on superconducting flux qubits have 
been reported 2 -* 3 -. However, these QND schemes work 
only far away from the degeneracy point (the 'sweet spot' 
where the sensitivity to noise is minimized), and the QND 
criterion is only approximately satisfied. Here we propose 
an ideal QND measurement scheme of a flux qubit that 
can be applied in a wide range of bias conditions. The 
QND fidelity for this measurement is expected to increase 
significantly as compared to previous proposals. 

In quantum mechanics, measurements induce back- 
action to the system under investigation due to Heisen- 
berg uncertainty. This back-action puts a fundamental 
limit on the precision of repeated measurements. In or- 
der to beat the standard quantum limits, the concept of 
QND measurement was developed in the context of grav- 
itational wave detection where repeated measurements 
beyond the standard quantum limit are required 5 . This 
concept has been extended from gravitational wave de- 
tection to other physical systems. A number of experi- 
ments has been performed in a micromechanical system^ 
and quantum optical systems 7 . This special type of mea- 
surement leaves the output state unaffected by subse- 
quent measurements and the free evolution of the system. 
QND measurements are crucial to overcome detector in- 
efficiencies, and to quantify the external disturbance to 
the QND variables. It is also found to have more versatile 
applications such as error correction 8 , one-way quantum 
computing 9 , low- noise amplification 10 and entanglement 
generatio n 11 ! 12 . 

In superconducting qubit systems, weak continuous 
QND measurements on superconducting transmon qubits 
have been realized in the dispersive limit 13 . This cir- 
cuit QED system has also been used to detect single mi- 
crowave photons in a coplanar wave guide 14 . Using the 
Josephson bifurcation amplifier, strong projective QND 
measurements have been demonstrated for quantronium 
qubits, flux qubits, and transmon qubits ^ 15 ! 16 . In order 
to implement QND detection for a continuous QND vari- 



able, a number of criteria have to be satisfied 4,5 . Among 
them, the most restrictive one is that the system free 
Hamiltonian H s commutes with the interaction Hi nt be- 
tween the system and the detector, i.e., [H Sl H int ] = 0. 
For existing flux qubit measurements 2,3 , this QND cri- 
terion is only approximately satisfied when the qubit is 
biased far away from the degeneracy point. However, 
the quantum coherence times for solid-state qubit vanish 
rapidly in this regime. QND detection close to the qubit 
degeneracy point is therefore desired. Moreover, to ac- 
quire full qubit control, the qubit bias has to be changed 
during various operations. After implementing an oper- 
ation at a certain bias, it is desirable to be able to carry 
out a QND measurement at that point, without adiabat- 
ically shifting back to another bias value. In this paper, 
by introducing an rf SQUID coupler to mediate the inter- 
action between a flux qubit and the detector, a Josephson 
bifurcation amplifier (JBA), we find a detection scheme 
that allows to implement a QND measurement at ar- 
bitrary bias including the degeneracy point. Moreover, 
our scheme works beyond the dispersive limit and can be 
extended to the case of strong qubit-detector coupling. 
This will help to improve the readout contrast to achieve 
a higher measurement fidelity and shorter measurement 
times. 

Another advantage of this scheme is the possibility to 
improve the so-called QND fidelity, which quantifies the 
accuracy of repeated measurements. In QND measure- 
ments by a JBA^ii 5 -, the drive on the JBA is first ramped 
to the bifurcation point to induce transitions between two 
bistable states. It is then reduced to maintain a latching 
plateau. The circuit geometry in the previous experi- 
ments does not implement an ideal QND measurement, 
i.e., [H Sl Hi nt ] 7^ 0. Qubit relaxation is then acceler- 
ated by the forced oscillations of the nonlinear resonator. 
The population fraction lost during the latching plateau 
and the preparation stage of the subsequent measure- 
ment limit the accuracy of the subsequent measurement. 
It turns out that the JBA induced qubit relaxation is 
the main limiting factor for the QND fidelity 17 . In our 
design, if proper control on the bias is acquired, the de- 
tection scheme is an ideal QND measurement. The QND 
fidelity is only limited by environment-induced qubit re- 
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FIG. 1: (Color online) Schematic diagram of the circuit. Red 
part: gradiometer flux qubit to be measured. Blue part: mea- 
surement device, a Josephson bifurcation amplifier formed by 
a dc SQUID shunted by a capacitance. Green part: rf SQUID 
acting as a tunable coupler. Orange parts: bias circuits. 



laxation, which is usually one order of magnitude smaller 
than the JBA induced relaxation. 

The structure of this paper is as follows: In Sec. II, 
we describe the circuit layout and the effective mutual 
inductance between the qubit and the JBA. The QND 
feature of this detection is analyzed in Sec. Ill, where two 
situations with bias at and off the degeneracy point are 
discussed respectively. In Sec. IV, we revisit the working 
principle of the JBA and calculate the qubit relaxation 
rates in the measurement process. With those relaxation 
rates, the fidelity of the QND measurement is evaluated. 
Section V discusses and summarizes our results. 



II. THE CIRCUIT LAYOUT 

Previous measurements can only work in the regime 
far away from the degeneracy point. This is because 
the measurement circuit (e.g. a Josephson bifurcation 
amplifier) can only be coupled to supercurrents in the 
loop. However in the conventional 3- Josephson junction 
desig n 18 ! 19 , the current states are the eigenstates of the 
system only if the qubit is biased to the degeneracy point. 
A natural solution for this problem is to use a gap-tunable 
qubit— - — and couple the measurement device to the dc 
SQUID part. This will enable a QND measurement when 
the qubit is biased at the degeneracy point. However, 
to implement a QND measurement at variable bias, the 
coupling with the measurement device has be to medi- 
ated in a way that it can always follow the eigenstates of 
the system. In this paper, we introduce a tunable coupler 
between the flux qubit and the Josephson bifurcation am- 
plifier. The qubit shares one control line with the tunable 
coupler. As the qubit bias is varied, the qubit coupling to 
the Josephson bifurcation amplifier is modified simulta- 



neously. We find that under certain conditions, a perfect 
QND measurement can be performed at variable qubit 
bias. 

The system we have in mind is shown is shown in 
Fig. [TJ It is composed of four parts: the system to 
be measured (red part), the measuring apparatus (blue 
part), the coupler (green part) and the bias circuits (or- 
ange parts). The system to be measured is a gradiometer- 
type superconducting flux qubi t 2Q i 21 which contains four 
Josephson junctions in three loops: The two lower loops 
(the main qubit loops) and the upper loop (the dc SQUID 
loop) penetrated by magnetic fluxes 3>i, $2, and $3. The 
two junctions in the dc SQUID loop are assumed to have 
identical Josephson energies a^Ej, here ao is the ratio 
between the Josephson energy of the smaller junctions 
and that of the two bigger junctions. The other two 
junctions are assumed to have the Josephson energy Ej. 
The total Josephson energy of the circuit is 2 ^ 

Ej cos (fi + Ej cos (p 2 + olEj cos (27r$ £ /$o - (<pi + ^2)) 

(1) 

where $ £ = ($1 - $2)/2, a = 2a cos(7r^ 3 /^> ), and 
(fk (fc = 1, 2, 3, 4) is the phase difference across the k-th 
Josephson junction. If <fr t is chosen close to 3>o/2 where 
$0 = h/(2e) is the flux quantum, the circuit dynam- 
ics can be effectively described in a two-level subspace 
of a double-well potential, and thus constitutes a flux 
qubit 18,19 . Together with the charging energy, the total 
Hamiltonian of the qubit is 



H = e($ t )<Tz + A($ 3 K 



(2) 



The Pauli matrices read a z = |0)(0| — cr x = 

1 0> < X I + |1)(0|, where |0) and |1) denote the states with 
clockwise and counterclockwise currents in the outer 
loop. The energy spacing of the two current states is 
s($> t ) = Ip($t — $o/2), with I p the magnitude of the clas- 
sical persistent current in the loop. The tunneling ampli- 
tude between the two states A ($3) = A (a) depends on 
the bias in the dc SQUID loop. Note that in contrast to 
the original flux qubit design 18 -^, this gradiometer flux 
qubit is insensitive to homogeneous fluctuations of the 
magnetic flux 20 . More importantly, it enables the JBA 
to couple with the dc- SQUID loop without changing the 
total bias flux of the qubit. 

The detector for the flux qubit is a Josephson bifur- 
cation amplifier (JBA) 24 (blue part in Fig. [1]), which in 
our scheme is a dc SQUID shunted by a capacitance C, 
subject to a microwave drive Irf cos (oodt + 4>a)- The 
JBA SQUID loop contains two identical Josephson junc- 
tions of critical current Iao- The phase differences across 
the two junctions are denoted by tpAi, ¥A2 respectively. 
The current in the loop is I a = Ia($a) cos</?a 5 with <&a, 
the flux bias in the JBA SQUID, set by external coils, 
Ia($a) = 2/ A osin(7r$A/$o), and <p A = (<p A1 + cp A 2)/2. 
The JBA circuit forms a driven nonlinear resonator which 
exhibits bistable behavior with hysteresis. With appro- 
priate drive sequences, a transition to one of the bistable 
states is correlated with the qubit states in a probabilis- 
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tic way. Therefore the qubit state can be read out by the 
phase of the transmitted or reflected microwave. 

The flux qubit is coupled to the JBA through their 
mutual inductance. There are two contributions to their 
mutual inductance: the direct mutual inductance (DMI) 
Majz (k = 1, 2, and 3 denotes the different loops in the 
qubit) and the effective mutual inductance (EMI) M' Ak . 
Hence the JBA produces flux biases to the qubit loops 
of the form (M Ak + M' Ak )I A . The EMI is induced by the 
nearby rf SQUID which acts as a coupler for the qubit 
and the JBA. The self-inductance of the coupler is as- 
sumed to be much larger than the mutual inductances 
and the dynamics of the coupler is confined to its lowest 
energy bands 25,26 . The DMI is fixed by fabrication pro- 
cesses while the EMI is tunable by the magnetic bias of 
the coupler <£>c as 2 ^ 2 ^ 



M' Ak ^c) 



M AC M C k P c cos0 



L 



c 



1 + (3 C cos #o 



(3) 



where Oq satisfies the nonlinear equation 



2tt 



= TiH^ + M AC I A o) ~ /3 c sin<9 (4) 



with f3 c = 2ttLcIco/&o, is the self-inductance of the 
coupler, and Ico the circulating coupler current. In par- 
ticular, if the coupler is biased at $>c = ((2n + 1/2)tt + 
f3 c )$o/27r-M AC I A o or ((2n+3/2)7r-f3 c )<$> /27r-M AC I A o 
(n is an arbitrary integer), the effective mutual induc- 
tance vanishes, M Ak = 0. Thus, for this bias condition, 
the EMIs between the JBA and all the qubit loops are 
canceled, and only the DMIs contribute to the coupling. 

Besides tunability, there is another important differ- 
ence between the DMI and the EMI: the DMI is sym- 
metric with respect to the qubit loops 1 and 2, while 
the EMI is not symmetrical, that is, M A i = M A 2, while 
M' Al ^ M' A2 (since M C i ^ M C2 ). Hence only the EMI 
couples the JBA to the gradiometer qubit flux & t in the 
form 



* t = (M' A1 - M' A2 )I A . 



(5) 



In our scheme, the whole chip is biased by external 
coils so that a homogeneous magnetic field threads all 
the loops. By choosing the area of each loop appropri- 
ately, the required background bias values can be im- 
posed. Besides the coupling to the external coils, the dc 
SQUID loop of the qubit is also coupled with an on-chip 
bias current Ib3 through Mb3- The qubit loop 1 shares 
another on-chip bias (the lower orange part) with the 
coupler. A bias current Ib in this bias line couples to the 
qubit loops and the coupler loop through mutual induc- 
tances Mbu and Mbc- In the following discussion, we 
will see that this shared bias is crucial for the possibility 
to do a QND measurement at arbitrary bias. 



III. QND NATURE OF THE DETECTION 
SCHEME 



In order to analyze the QND nature of the detection 
scheme, we first look at the situation that the bias is 
set at the degeneracy point and then study the case of a 
general (off-degeneracy) bias. 



A. Degeneracy point 

We first look at the case when the qubit is biased at 
the degeneracy point <& t b = &o/2. At this point, the first- 
order flux noise disappears so that the qubit quantum 
coherence can be preserved longer. 

The qubit is biased at the degeneracy point by trapping 
one flux quantu m 20 ! 21 . The bias current is set to zero 
Ib = and the flux bias of the coupler is set by external 
coils to be 



Cb 



$0 

2tt 



(§+*■) 



M AC I A ■ 



(6) 



According to the discussion following Eq. (j3j), the effec- 
tive mutual inductance M Ak vanishes at this bias. Thus, 
the qubit only couples to the JBA through the direct 
mutual inductance. As shown in Eq. (j5]), this means the 
JBA has no influence on but only couples to $3. If 
7tM A3 I a <C $0? the Hamiltonian can be expanded to first 
order as 23 



Ha — H ( 



qO 



(7) 



where H q o is the free Hamiltonian of the qubit 

H q0 = A($3bK , (8) 
and Hj is the interaction between the qubit and the JBA 
Hi = \($3b)crx cos (p A . (9) 

where $35 is the total flux bias of the dc SQUID loop 
(generated by both external coils and I 33)- The coupling 
coefficient is 



A($36) 



ttM A3 I a 

$0 



K($3b) 



(10) 



with 



K(®3b) = 2a sin ( tt 



$36 \ dA(a) 



da 



and a = 2ao cos^^b/^o). The coupling energy be- 
tween the qubit and the JBA is tunable by $35. 

Equations (0 - (j9j) show that the free Hamiltonian 
commutes with the interaction Hamiltonian. If one 
chooses a x as the measurement observable, a QND mea- 
surement can be implemented. 
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B. General (off-degeneracy) bias 

If we change the current in the shared bias by a small 
amount 1b = the qubit is biased away from the 

degeneracy point, and the corresponding bias change in 
the qubit loop is 5$> t = (M B i — Mb2)5I B - Since the 
coupler shares the same bias, the magnetic flux pene- 
trating the coupler bias is also shifted by a small amount 
MbcSIb (MbcSIb <C &cb is always satisfied in the rel- 
evant operation regime). This will induce a non-zero 
effective mutual inductance M' Ak (<&cb + Mbc$1b) = 
-(27r/$o) 2 I c oMAcMckMBc5lB. As we discussed after 
Eq. (j5j), a non-zero EMI will couple the JBA to the qubit 
flux bias & t as weu as $3- The qubit Hamiltonian under 
this bias reads 



with 



H q — H q o + Hi 



H q o = Vt z (J z + £l x Vx 



and the interaction Hamiltonian 

Hi = (X z cr z + \ x ($3b)crx)cos(pA 
with Q z = I p 8$ t , Q x = A($ 3b ) and 

X z = I P (M' A1 - M' A2 )Ia , 



(11) 
(12) 

(13) 

(14) 
(15) 



If we define a parameter 77 = \ z Vt x / \ x Vt Zl it is straight- 
forward to see that the free qubit Hamiltonian commutes 
with the interaction Hamiltonian when 77 = 1. Therefore, 
a sufficient condition to implement a QND measurement 
at variable flux bias is 



MbcIco Ma M AC A($ 3b ) 
$ M B i Mas ^sb) 



= 1 



(16) 



where we have neglected M#2, Mc2, and Mc3 since they 
are much smaller than the other mutual inductances. 

In other words, if 77 = 1, as the qubit is biased away 
from the degeneracy point, the interaction with the JBA 
is changed accordingly, so that the interaction Hamilto- 
nian always commutes with the qubit free Hamiltonian. 
This condition is possible to be satisfied experimentally, 
e.g., if ao = 0.4, and the bias $35 satisfies a = 0.7. At 
this bias, A/(dA/da) « -0.11 28 . If I C q = 1 /xA, M B c = 
23.5 pH, Mci = 25 pH, M B i = 5 pH, Mac = 25 pH, 
Ma3 — 5 pH, then 77 = 1. Note that 77 depends on the 
bias $35 which is tunable in situ, i.e., by tuning Ib3- 
Thus, errors in the fabrication process can be compen- 
sated to satisfy Eq. ([16]) . the condition for QND detec- 
tion. Note that $35 is determined by this condition since 
all the other parameters are fixed by fabrication. There- 
fore, this QND scheme works for variable bias values of 



the main qubit loop, but it does not work for variable bias 
values of the dc SQUID loop in the general case. How- 
ever, we would also like to point out that if the qubit is 
biased away from the degeneracy point, by changing e, 
any arbitrary single-qubit operation can be implemented; 
in this sense, it is not necessary to tune $35. Also, if the 
qubit is biased at the degeneracy point, which is also the 
situation in which tuning $35 is meaningful, the QND 
measurement can be implemented for variable $35. 

For 77 = 1, the results obtained for the degeneracy 
point and general (off-degeneracy) bias can be written in 
a uniform way as 



H ( 



qO 



Qcr z 



with ft = + Of, and a z = (Q z cr z 

interaction Hamiltonian is 

Hi — Xa z cos if a 



(17) 

n x a x )/n. The 
(18) 



with A = \Jy? z + X x . At the degeneracy point, Vt z = X z = 
0, so that a z = <j x . 

IV. MEASUREMENT FIDELITY 

The JBA is an oscillator with nonlinear Josephson in- 
ductance. Under a strong microwave drive, the Joseph- 
son energy of the junction —Eja cosipA is expanded be- 
yond the harmonic approximation and the classical dy- 
namics can be described by a Duffing oscillator—. For a 
certain range of drive conditions, the nonlinear oscillator 
exhibits bistable behavior with hysteresi o 24 i 29 . The two 
possible stable states correspond to different oscillation 
amplitudes and phases, which can be distinguished by 
transmitted or reflected microwave signals 2 -^^ . Switch- 
ing between the two stable states happens when the 
drive power reaches a certain threshold. The switch- 
ing probability depends on the value of the nonlinear 
inductance, which in our case depends on the states of 
the qubit through the mutual inductance. This is be- 
cause the effective Josephson energy of the junctions 
of the JBA is modified by the interaction Eq. ([15)) as 
Eja(&z) = -?A$o/27r — \(J Z - Therefore, measuring the 
phase of the transmitted microwave signal, one can read 
out the state of the qubit. 

The back-action from the measurement device destroys 
the phase coherence of the qubit states during the read- 
out process. Besides dephasing, the back-action could 
also induce relaxation to the qubit. This is the case for a 
qubit Hamiltonian H q o = £lcr z + Aa x with a small non- 
ideal QND fraction A, see e.g. the QND measurement of 
Ref. where A/Q w 0.34. The JBA is strongly coupled 
to a dissipative environment while weakly coupled to the 
qubit. Hence it serves as a bath for the qubit. According 
to Eq. ([15]) , the influence of the JBA on the qubit can be 
described by its correlation function 



f 

Jo 



dte ZUJt (cos (pA(t) cos(/?a(0)) , 



(19) 
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and the induced decay rate can be calculated through the 
Fermi golden rule. The Bloch Redfield rates induced by 
the operation of the JBA are 



r M 



A 2 A 



^9? + A 2 

\ 2 n 



&(G(y/n 2 + A 2 )) 



Vfl 2 + A 2 



5=»(G(0)), 



(20) 



where T r is the induced relaxation rate and is the in- 
duced dephasing rate. When the JBA is ramped to the 
measurement level and the latching plateau, the correla- 
tion function is prominently increased due to quantum 
activation^. Qubit decay is enhanced by the measure- 
ment operatio n 17 ! 31 . This results in qubit relaxation and 
the measurement is driven away from the QND regime. 
This induced relaxation has been found to be the main 
source of measurement error^i 7 -. One way to reduce this 
back-action is working in the dispersive limh— . In our 
case, an ideal QND measurement is possible, i.e. A = 0, 
so that T r = 0, i.e., the JBA does not induce extra relax- 
ation but only dephasing to the qubit. Hence the QND 
condition can be preserved better in our scheme and the 
QND fidelity can be improved. 

Besides the induced decay rates Eq. (|20|) . there is an- 
other decay mechanism due to the flux fluctuations of the 
environment. This will perturb the fluxes in the qubit 
loops as 5& t — fJ>tX and S^s = [I3X, where X represents 
an environmental operator (such as a two-level fluctua- 
te) and fi t (fjis) characterizes its coupling strength to the 
different qubit loops. Hence the qubit is coupled to the 
environment as AH = ^tXa z + £3X0^, with £ t = I p ^op t 
and £3 = — 7r/i3^(^>36)/^>o- In the interaction picture 
AHj = (£ 3 sinx + £tCos X )X(t)4 ( 2 i) 
+ (6cos X -6sin X W) (<j;e* m + <7'_e-* m ) , 

with cosx = e/Cl and sinx = A/O. 

According to the Fermi golden rule, the relaxation 
rate is = (£3 cosx — ^ sm x) 2 ^x(^ = ±0), where 
S x (w) = J^dr{X(r)X({)))e iuJT is the flux noise spec- 
trum. In a real experiment, the flux noise could have 
multiple sources, such as two- level fluctuators inside the 
barrier, high-frequency noise from the control lines 20 , and 
others. Therefore the noise spectrum may exhibit a com- 
plicated frequency dependence and have a strong sam- 
ple dependence. In our discussion, we assume an Ohmic 
noise spectrum (/-noise) for the environment bath plus 
a few peaks due to two-level fluctuators 

Sx(ou) = cuR (^coth (^j^j^j + 1 ) + ^2S i 5(uj-uj i ), 

(22) 

where Rq is the Ohmic impedance and To is the environ- 
mental temperature. The QND fidelity of two successive 
measurements is 



^QNd(t) 



p(e\e) +p(g\g) = exp(-iy) + exp(~r t r) 
2 2 

(23) 



F — 




FIG. 2: (Color online) Dependence of the QND fidelity de- 
fined in Eq. (|23|) on the bias of the flux qubit for different 
values of k = Inset: dependence of the QND fidelity 

on the bias flux (in units of &o/tt) in the dc SQUID loop, 
when the qubit is biased at the degeneracy point. 



where p(e\e) (p(g\g)) is the probability that the qubit 
state \e) (\g)) is unchanged after the first measurement 
and t is the time interval between the two measurements. 

Figure [2] shows the dependence of the QND fidelity on 
the qubit bias e for different values of k = £3/^ (usu- 
ally k < 1 since the perturbation on the main qubit loop 
in general has a larger influence than the perturbation 
on the dc SQUID loop 28 ). Here we assume that the 
time interval between two measurements is r = 50 ns, 
To = 20 mK, the qubit relaxation time is 250 ns at the 
degeneracy point, and the value of A is fixed at 7.8 GHz. 
The plot shows that the measurement fidelity remains 
rather high for a wide range of bias values. Even at the 
degeneracy point e = where the relaxation is strong, 
a measurement fidelity larger than 90% can be achieved. 
As the bias is increased to the positive side, the fidelity 
increases as the relaxation decreases. Far above the de- 
generacy point, the measurement fidelity is very close to 
100%. Note that the fidelity is not symmetrical with re- 
spect to the axis e = but becomes more symmetrical 
as k decreases. At k = 0, the curve shows complete sym- 
metry because noise only contributes to the main qubit 
loop, it is symmetrical with respect of the sign of the 
qubit bias. The inset of Fig. [2] shows that the QND fi- 
delity at the degeneracy point decreases with the SQUID 
bias Tr^sb/^o- This is because at the degeneracy point, 
the qubit relaxation rate due to /-noise increases linearly 
with the gap A and A increases with the SQUID bias 28 . 
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FIG. 3: (Color online) Schematic diagram of the circuit for a 
non-gradiometer qubit. Red part: non-gradiometer flux qubit 
to be measured. Blue part: measurement device, a Josephson 
bifurcation amplifier formed by a dc SQUID shunted by a ca- 
pacitance. Green part: rf SQUID acting as a tunable coupler. 
Orange parts: bias circuits. 



FIG. 4: (Color online) Dependence of the scaled effective mu- 
tual inductance M' — M' Ak / '{MacMch I ' Lc) on the coupler 
bias <I> = &cb+MacIao, see Eq. J3]). The red dots on the two 
sides (intersections with the dashed line) indicate the back- 
ground bias in the case of the gradiometer qubit, while the 
green point in the middle is an example of the background 
bias for a non-gradiometer qubit. Here /3 C = 0.9. 



The measurement fidelity can be used as a noise spec- 
trometer for environmental fluctuations. This is actually 
one of the main applications of QND measurements: de- 
tecting perturbations to the system. The QND nature 
of the measurement guarantees that the readout back- 
action will not change the value of the observable. The 
measurement fidelity therefore reflects the noise spec- 
trum of the environment. For example, the existence 
of one two-level fluctuator inside the barrier 32,33 would 
be revealed by a corresponding peak in the QND fidelity 
at a certain bias. 



V. DISCUSSION AND SUMMARY 

All the discussion above is based on the gradiometer 
type- flux qubit. However, with a few modifications as ex- 
plained below, the measurement protocol can be adapted 
to non-gradiometer type flux qubits with a tunable gap 22 , 
see Fig. [3j Two current bias lines Icb and lib are used 
to control the coupler and the qubit separately in order 
to guarantee a QND measurement for a non-gradiometer 
flux qubit at the degeneracy point. 

The background bias of a non-gradiometer qubit is sen- 
sitive to homogeneous magnetic field fluctuations, but 
can be implemented easily by external coils (while the 
gradiometer qubit requires the technique to trap flux- 
oids). Also, it is possible to achieve a more sensitive tun- 
ing comparing with the gradiometer qubit. This can be 
seen from Fig. 2] which shows the scaled effective mutual 
inductance with respect to the coupler bias. To achieve a 
more sensitive tuning within the tunable range of the on- 
chip bias (typically on the order of m$o = 10 _3< l ) o), the 
coupler is desired to be pre-biased close to 3>o/2. How- 



ever, in the case of the gradiometer qubit, since M' Ak 
should be zero when the qubit is biased at the degener- 
acy point, the background bias should be set around the 
red points in Fig. [4j i.e., relatively far from $>o/2. For 
non-gradiometer qubits, the background bias point is de- 
termined by the fabrication process. If the mutual induc- 
tance between the JBA and qubit is large and the coupler 
loop is twisted as indicated in Fig.[4j the background bias 
can be set closer to <&q/2 (say, the green point in Fig. 2]). 
As a result, a more sensitive tuning can be achieved. For 
example, if we assume Ma3 = 5 pH, Mai — 0.5 pH, 
M A c = 10 pH, L c = 100 pH, M C i = 10 pH, I A = 1 /iA, 
for a change in the qubit bias 5$> t = 2 m$o, the coupler 
bias &cb should be tuned by 15 m^o in the case of the 
gradiometer qubit, while only 0.4 m$o in the case of a 
non-gradiometer qubit. 

In summary, we have studied a scheme to real- 
ize a quantum non-demolition (QND) measurement for 
gradiometer-type flux qubits by a Josephson bifurcation 
amplifier. We have shown that a perfect QND mea- 
surement can be implemented for a qubit with variable 
magnetic bias. The QND fidelity of this measurement 
is expected to be high over a wide range of bias con- 
ditions. We have also discussed how to generalize our 
scheme to non-gradiometer qubits. Our estimates indi- 
cate that such a QND measurement may be realized ex- 
perimentally, and we hope that this will happen in the 
close future. 
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